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A randomly interacting N-species Lotka-Volterra system in the pres- 
ence of a Gaussian multiplicative noise is analyzed. The investigation is 
focused on the role of this external noise into the statistical properties of 
the extinction times of the populations. The distributions show a Gaussian 
shape for each noise intensity value investigated. A monotonic behavior of 
the mean extinction time as a function of the noise intensity is found, while 
a nonmonotonic behavior of the width of the extinction time probability 
distribution characterizes the dynamical evolution. 

PACS numbers: 05.40.-a, 05.45.-a, 87.23.Cc, 89.75.-k 

1. Introduction 

Generalized Lotka-Volterra equations have been used in recent years 
to describe the dynamics of various kind of population species, which are 
main components in complex ecosystems ^J-^- To understand the com- 
plex behavior of such ecosystems is crucial to analyze the role played by 
the external noise on the dynamics. It has become increasingly evident 
that nonlinearity and noise play an important role in such complex dy- 
namics. Recently, in fact, noise-induced effects in population dynamics, 
such as pattern formation, stochastic resonance, noise delayed extinction, 
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quasi periodic oscillations etc... have been investigated with increasing in- 
terest [8j-[14J. Complex ecological systems evolve towards the equilibrium 
states through the slow process of nonlinear relaxation, which is strongly 
dependent on the random interaction between the species, the initial con- 
ditions and the random interaction with environment. One of the open 
problems of such ecosystems is the investigation of the time scales of extinc- 
tion and survival of the species and their related statistics. Various factors 
affecting extinction such as migration, chaos, interaction between species, 
spatial synchronization, etc., have been discussed in the literature [15j-|18j. 
However, there is lack of investigation on the role of external noise on the 
extinction process, which is the main focus of this paper. The mathematical 
model used to analyze the dynamics of biological species, with spatially 
homogeneous densities, is the generalized Lotka-Volterra system. We con- 
sider a Malthus-Verhulst modelling of the self regulation mechanism and an 
external multiplicative noise source, taking the environment interaction into 
account \19\ [20] . Within this model we analyzed the role of the noise in the 
statistical properties of the extinction times of the populations. Specifically 
a monotonic behavior of the mean extinction time as a function of the noise 
intensity is observed. The width of the distribution of the extinction times, 
however, has a nonmonotonic behavior as a function of the noise intensity. 



2. The model 

The dynamical evolution of our ecosystem composed by N interacting 
species in a noisy environment (climate, disease, etc..) is described by the 
following generalized Lotka-Volterra equations with a multiplicative noise, 
in the framework of Ito stochastic calculus |21j 



dni{t) 



!a + - ) -ni{t) + 



Jijnj{t) I dt + ^/edwi 

I 



riiit), (1) 



where nj(t) > is the population density of the i*^ species at time t and 
i = 1, ...,N. In Eq. ([1]), the first two terms describe the development of the 
i^^ species without interacting with other species, a is the growth parameter, 
and Jij is the interaction matrix, which models the interaction between 
different species {i ^ j). Here Wi is the Wiener process whose increment dwi 
satisfy the usual statistical properties {d'Wi{t)) = 0, and {dwi{t)dwj{t')) = 
8ij5{t — t')dt. The interaction matrix Jij has elements randomly distributed 
according to a Gaussian distribution with {Jij) = 0, {JijJji) = 0, and a'j = 



J /N. Our ecosystem contains, therefore, 50% of prey-predator interactions 
{Jij < and Jji > 0), 25% competitive interactions {Jij > and Jji > 0), 
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and 25% symbiotic interactions {Jij < and Jji < 0). We consider all 
species equivalent so that the characteristic parameters of the ecosystem 
are independent of the species. The formal solution of Eq. ([1]) is 



ni{0)exp 


at + ^/ewi{t) + /q dt' Y.jjti Jij^jit')) 




1 + nj(0) /q dt'exp 


at' + ^/ewi{t') + /q dt" Y.j^i Jij"^ 





(2) 



where the term hi(t) = X^j^j </^j?^j(^) represents the influence of other 
species on the differential growth rate of the i^^ population. The dynami- 
cal behavior of the i*^ population depends on the time integral of the term 
hi(t) and the time integral in the denominator of Eq. ([2|). By considering 
the deterministic dynamics (in the absence of external noise (e = 0)), with 
a large number of interacting species (that is large interaction random ma- 
trix), we can assume that the term hi(t) is Gaussian with zero mean and 

variance cr^^ = llj^k{JijJik){njnk) = J'^{nf), with {JijJik) = Sjkir- In the 
absence of external noise, from the fixed-point equation nj(a — nj + hi) = 0, 
the stationary probability distribution of the populations is the sum of a 
truncated Gaussian distribution at nj = (nj > 0) and a delta function for 
extinct species. The initial values of the populations nj(0) have also Gaus- 
sian distribution with mean value (rej(O)) = 1, and variance c^(o) ~ 0.03. 
The interaction strength J determines two different dynamical behaviors 
of the ecosystem. Above a critical value Jc = 1-1, the system is unstable 
and at least one of the populations diverges. Below Jc the system is stable 
and asymptotically reaches an equilibrium state. The equilibrium values of 
the populations depend both on their initial values and on the interaction 
matrix. If we consider a quenched random interaction matrix, the ecosys- 
tem has a great number of equilibrium configurations, each one with its 
attraction basin. For an interaction strength J = 1 and an intrinsic growth 
parameter a = 1 we obtain: (n^) = 1.4387, (n?) = 4.514, and cr^. = 2.44. 
These values agree with that obtained from numerical simulation of Eq. ([1]). 

The statistics of the species extinction has been analyzed using the mean 
extinction time (tm), defined as 



^ ^ exp 

{tm) = -TJ tm, (3) 

and its variance 



= {tl) - {tm)'. 



(4) 
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Here {tm) is an ensemble average {N^xp is the number of simulative exper- 
iments), tm is the average extinction time over the number of populations 



i=l 

and ti^m is the extinction of the i-th population in the m-th experiment. 

3. Results and Comments 

The parameters used in our simulation are: a = 1.2, J = 1, ctj = 
0.005, = 400. The number of simulative experiments is N^^p = 1000, 
and the initial values of the average population and its standard deviation 
are: (jij) = 1, cr^^ = 0.03. The dynamics of various species are different 
even if they are equivalent according to the parameters in the dynamical 
equation ([1]). However to change the species index by fixing the random 
matrix or to change the random matrix by fixing the species index should 
be equivalent as regards the asymptotic dynamical regime. 

Probability of the Density Species: P(n j) 
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Fig. 1. Long time probability distribution of the species densities for different 
external noise intensities e. Namely e = 0, 0.562, 1, f. 778, 3.162, 5.62. Around the 
value e — 0.562 the distribution becomes asymmetric, and for e > 1.778, all the 
species are extinct. 
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In the presence of external noise (e 7^ 0) we calculate the long time 
probability distribution of the species density for different values of the 
noise intensity. These are shown in Fig. [TJ 

For increasing external noise intensity we obtain a larger probability dis- 
tribution with a lower maximum (see the different scales in Figs. [J for dif- 
ferent noise intensity values). The distribution is asymmetric for e = 0.562 
and tends to become a truncated delta function around the zero value 
(P(nj) = 5{ni) for rij > 0, and P{ni) = for rii < 0), for further increas- 
ing noise intensity. Specifically for high values of noise intensity (namely 
for e > 1.778) we strongly perturb the population dynamics and because 
of the presence of an absorbing barrier at rij = [T9], we obtain quickly 
the extinction of all the species. To confirm this picture we calculate the 
time evolution of the average number of extinct species for different noise 
intensities. This time behavior is shown in Fig. [2l We see that this number 
increases with noise intensity, obtaining a rapid transient dynamics of the 
system towards the extinction final state for e > 1.778. This means that the 
species rapidly die and the probability distribution of the species density 
confines accordingly. 

Ratio of Extinct Species: N f(t)/N 
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Fig. 2. Time evolution of the normalized number of the extinct species for different 
noise intensities e. Namely: e = 0, 0.562, 1, 1.778. 

In the following Fig. [3] we show the probability distribution function 
(PDF) of the extinction times of the species. The shape of the distribution 
is Gaussian in the deterministic regime (e = 0) and in the presence of the 
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external noise (e 7^ 0). For low noise intensities the probability distribu- 
tion becomes larger and lower until reaches the value of e = 1. After this 
value of noise intensity the distribution becomes narrow and higher. The 
mean extinction time, which is easily visible from the figure because of the 
Gaussian shape distribution, decreases monotonically with increasing noise 
intensity. In this figure it is shown a well defined extinction time windows of 
the species, moving towards the absorbing barrier at rii = 0, with increasing 
noise intensity. 
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Fig. 3. Probability distribution function of the extinction times of the species, for 
different values of the noise intensity. Namely: e = 0,0.562,1,1.778,3.162,5.62. 
All the PDFs show a Gaussian shape distribution. 

This behavior is due to the presence of the external noise and the ab- 
sorbing barrier. In fact, in the deterministic case (e = 0), the Gaussian 
distribution of the extinction times is due only to the random interaction 
matrix. The characteristic values of the distribution, that is the mean and 
the variance, depend on the choice of the parameters of the model, that is 
the growth parameter a, the interaction strength J and the initial condi- 
tions. A small amount of noise forces the system to sample more of the 
available range in the parameter space and therefore moves lightly the sys- 
tem towards the extinction. The average extinction time at e = 0.562 is 
less than that in the absence of external noise. This enlargement and low- 
ering of the PDF continues until the noise intensity reaches the value of the 
interaction strength J = 1. After that the external noise prevails on the in- 
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teraction matrix term and the extinction process proceeds quickly because 
of the presence of the absorbing barrier at rij = (see Eq. ([I])). Almost 
all the species extinguish in short times around a very low mean extinction 
time. At e = 3.162, for example, (tm) ~ 3.3. Increasing the noise intensity 
(e > 1), therefore, the PDF becomes narrower and higher. 

As can be seen in Fig. [3]for e = 3.162 and e = 5.62 the species extinction 
happens in few time units, so that the probability of density species vanishes 
for the same values (see Fig. [T]). 

This peculiar behavior of the PDF of extinction times gives rise to the 
nonmonotonic behavior of the variance of the same quantity as a function 
of the noise intensity. This is shown in the following Fig. HI 
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Fig. 4. Mean extinction time and variance as a function of the noise intensity e. 
The variance shows a nonmonotonic behavior with a maximum at e ~ 1 and very 
low values at higher noise intensities. 
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From this figure, the maximum of the variance at the noise intensity 
e ~ 1 and the very small values of the variance at high noise intensities, 
are clearly visible. In the same figure the monotonic behavior of the mean 
extinction time {tm) is shown. Calculation have been repeated for different 
number of populations, namely: N = 100, 200, 300, 400. In all the calcula- 
tions the qualitative behaviors of the mean extinction time and its variance 
are the same than those reported in FigSl 

We didn't reveal any power law decay for the probability distribution 
of extinction times as found in previous investigations [16], but we plan 
to make a more detailed investigation of species lifetime distribution in a 
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forthcoming paper. 



4. Conclusions 

The analysis of the dynamics of ecosystem composed by N random in- 
teracting species has been performed in the presence of multiphcative noise. 
The probabihty density of the extinction time of the species {P{t)) has been 
evaluated for various noise intensities. The extinction times t^s are Gaus- 
sian distributed with a mean value monotonically decreasing as a function 
of the noise intensity. The variance of the extinction times shows a non- 
monotonic behavior, which characterizes the transient dynamics of the N 
random interacting species model. 
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